9 Quantization of Gauge Fields 



We will now turn to the problem of the quantization of the simplest gauge theory, 
the free electromagnetic field. This is an abelian gauge theory. In Physics 583 
we will discuss at length the quantization of non- abelian gauge fields. Unlike 
abelian theories, such as the free electromagnetic field, even in the absence of 
matter fields, non-abelian gauge theories are not free fields and have highly 
non-trivial dynamics. In Physics 582 we will discuss canonical quantization and 
path-integral quantization of the free electromagnetic gauge field. 



9.1 Canonical Quantization of the Free Electromagnetic 
Field 

Maxwell's theory was the first field theory to be quantized. However, the quan- 
tization procedure involves a number of subtleties not shared by the other prob- 
lems that we have considered so far. The issue is the fact that this theory has 
a local gauge invariance. Unlike systems which only have global symmetries, 
not all the classical configurations of vector potentials represent physically dis- 
tinct states. It could be argued that one should abandon the picture based on 
the vector potential and go back to a picture based on electric and magnetic 
fields instead. However, there is no local Lagrangian that can describe the time 
evolution of the system now. Furthermore is not clear which fields, E or B (or 
some other field) plays the role of coordinates and which can play the role of 
momenta. For that reason, one sticks to the Lagrangian formulation with the 
vector potential as its independent coordinate-like variable. 
The Lagrangian for Maxwell's theory 

C = -\F^ (1) 

where F^ v = d^A u — d v A^, can be written in the form 

C= 1 -(E 2 -B 2 ) (2) 

where 

Ej = -OqAj — djAo 
Bj = -ejkedkAt 

(3) 

The electric field Ej and the space components of the vector potential Aj form 
a canonical pair since, by definition, the momentum H, conjugate to Aj is 



dC 

SdoAj (x) 



n i( X ) = ra a — = 9 A 3 + d 3 A = ~ E 3 ( 4 ) 
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Notice that since L does not contain any terms which include <9o^4o, the mo- 
mentum IIo, conjugate to Aq, vanishes 

n = j^t- = o (5) 

5d A 

A consequence of this result is that A is essentially arbitrary and it plays 
the role of a Lagrange multiplier. Indeed it is always possible to find a gauge 
transformation <j> 

A' = A + d <f> A'^Aj-drf (6) 
such that A' = 0. The solution is 

d <j> = -Ac (7) 

which is consistent provided that Aq vanishes both in the remote part and in 
the remote future (t — > ±oo). 

The canonical formalism can be applied to Maxwell's electrodynamics if we 
notice that the fields Aj(x) and Uj'(x') obey the equal-time Poisson Brackets 

{A j (x),n j ,(x*)} PB =6 jj ,6 3 (x-x') (8) 

or, in terms of the electric field E, 

{A j (x),E f (x')} PB = -S jr 5 3 (x-x') (9) 

The classical Hamiltonian density is defined in the usual manner 

H = n j d A j -£ (10) 

We find 

H(x) = ±(E 2 + B 2 )-A (x)v-E(x) (11) 

Except for the last term, this is the usual answer. It is easy to see that the last 
term is a constant of motion. Indeed the equal-time Poisson Bracket between 
the Hamiltonian density H(x) and V • E(y) is zero. By explicit calculation, we 
get 

{H{z),s?-E(3,)}pB-Jdz[- SMg) 5Ej{g) + 5Ej{g) SAj{zl ] (12) 
But 

jAM = J d W lBm SAM = Jd 3 wB k ( W )6(x- W )e kej v7S( W -z) 

= -ekljVe J d 3 wB k (w)5(x - w)S(w - z) 

(13) 
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Hence 

5H(x) 



5A 3 {z) 
Similarly, we get 



ejek Vl {B k (x)S(x- z)) = e jik B k (x) y? 5(^-5) (14) 



Thus, the Poisson Bracket is 



V ' ^(y)}ps = y d 3 z[-e jek B k (x) V ? 5(2-2) V V j S(y- z)] 

= -e jtk B k (x) y? V V j5{x-y) 
= e jek B k (x) v? S7jS(x- y) = 

(16) 

provided that B(x) is non-singular. Thus, y • £?(a?) is a constant of motion. It 
is easy to check that y • i? generates infinitesimal gauge transformations. We 
will prove this statement directly in the quantum theory. 

Since y • E(x) is a constant of motion, if we pick a value for it at some initial 
time x = t , it will remain constant in time. Thus we can write 

V • E(x) = p(x) (17) 

which we recognize to be Gauss's Law. Naturally, an external charge distribu- 
tion may be explicitly time dependent and then 

l(V^) = |(V^) = |pe^,*) (18) 

Before turning to the quantization of this theory, we must notice that Aq plays 
the role of a Lagrange multiplier field whose variation forces Gauss's Law, y -E = 
0. Hence Gauss's Law should be regarded as a constraint rather than an equation 
of motion. This issue becomes very important in the quantum theory. Indeed, 
without the constraint y • E = 0, the theory is absolutely trivial, and wrong. 

Constraints impose very severe restrictions on the allowed states of a quan- 
tum theory. Consider for instance a particle of mass m moving freely in three 
dimensional space. Its stationary states have wave functions ^p{r, t) 



i 



p- f — E(p)t 



e V n J (19) 



p 2 

with an energy E(p) = - — . If we constrain the particle to move only on the 
2m 

surface of a sphere of radius R, it becomes equivalent to a rigid rotor of moment 
of inertia I = mR 2 and energy eigenvalues ee m = jj£(£+1) where € = 0,1,2,..., 
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and |m| < I. Thus, even the simple constraint r 2 = i? 2 , does have non-trivial 
effects. 

The constraints that we have to impose when quantizing Maxwell's electro- 
dynamics do not change the energy spectrum. This is so because we can reduce 
the number of degrees of freedom to be quantized by taking advantage of the 
gauge invariance of the classical theory. This procedure is called gauge fixing. 
For example, the classical equation of motion 

OA** - d ft (d v A v ) = (20) 
becomes, in the Coulomb gauge A = and y • A = 0, 

UAj = (21) 
However the Coulomb gauge is not compatible with the Poisson Bracket 

{Aj{S),TLj, (x')} PB = 6 jf 5(x - a?) (22) 

since the spatial divergence of the <5-function does not vanish. It will follow that 
the quantization of the theory in the Coulomb gauge is achieved at the price of a 
modification of the commutation relations. Because the classical theory is gauge 
invariant we can always fix the gauge without any loss of physical content. The 
procedure of gauge fixing has the attractive that the number of independent 
variables is greatly reduced. A standard approach to the quantization of a 
gauge theory is to fix the gauge first, at the classical level, and to quantize later. 
However, a number of problems arise immediately. For instance, in most gauges 
symmetries such as Lorentz invariance are lost, at least manifestly so. Thus, the 
Coulomb gauge, also known as the radiation or transverse gauge, spoils Lorentz 
invariance, but it has the attractive that the nature of the physical states (the 
photons) is quite transparent. We will see below that the quantization of the 
theory in this gauge has some peculiarities. 

Another standard choice is the Lorentz gauge 

d^A* = (23) 

whose main appeal is its manifest covariance. The quantization of the system is 
this gauge follows the method developed by and Gupta and Bleuer. While highly 
successful, it requires the introduction of states with negative norm (known 
as ghosts) which cancel-out all the gauge-dependent contributions to physical 
quantities. 

More general covariant gauges can also be defined. A general approach 
consists not on imposing a rigid restriction on the degrees of freedom, but to 
add new terms to the Lagrangian which eliminate the gauge freedom. For 
instance, the modified Lagrangian 

C = ~F a + ^{d^{x)f (24) 

is not gauge invariant because of the presence of the last term. We can easily 
see that this term weighs gauge equivalent configurations differently and the a 
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plays the role of a Lagrange multiplier In fact, in the limit a — > oo we recover 
the Lorcntz gauge condition. In the path integral quantization of Maxwell's 
theory it is proven that this approach is equivalent to an average over gauges of 
the physical quantities. If a = 1, the equations of motion become very simple, 
i.e., d 2 a M = 0. This is the Feynman gauge. This is the gauge in which the 
calculations are simplest. Still within the Hamiltonian or canonical quantization 
method, a third approach has been developed. In this approach one fixes the 
gauge A = 0. This condition is not enough to eliminate the gauge freedom. 
In this gauge a residual set of gauge transformations are still allowed, the time- 
independent ones. In this approach quantization is achieved by replacing the 
Poisson Brackets by commutators and Gauss' Law condition becomes now a 
constraint on the space of physical quantum states. So, we quantize first and 
constrain later. 

In general, it is a non-trivial task to prove that all the different quantizations 
yield a theory with the same physical properties. In practice what one has to 
prove is that these different gauge choices yield theories whose states differ 
from each other at most by a unitary transformation. Otherwise, the quantized 
theories would be physically inequivalent. In addition, the recovery of Lorentz 
invariance may be a bit tedious in some cases. There is however, an alternative, 
complementary, approach to the quantum theory in which most of these issues 
become very transparent. This is the path-integral approach. This method 
has the advantage that all the symmetries are taken care of from the out set. 
In addition, the canonical methods encounter very serious difficulties in the 
treatment of the non-abelian generalizations of Maxwell's electrodynamics. 

We will consider here two canonical approaches: 1) quantization in the 
Coulomb gauge and 2) canonical quantization in the Aq = gauge in the 
Schr'odinger picture. 

9.2 Coulomb Gauge 

Quantization in the Coulomb gauge follows the methods developed for the scalar 
field very closely. Indeed the classical constraints A = and V • A = allow 
for a Fourier expansion of the vector potential A(x,x ). In Fourier space we 
write 

f d 3 p 

A(x, x ) = J ^3 2 ^ A(p, x ) cxp(ip- x) (25) 

where A(p, xo) = A*(—p,x ). Maxwell's equations yield the classical equation 
of motion 

\ji(S,x )=0 (26) 
The Fourier expansion is consistent only if A(p, xo) satisfies 

d 2 A(P,x )+p 2 A(p,x ) = (27) 
The constraint y • A = now becomes the transversahty condition 

p-A(p,x )=0 (28) 
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Thus, A(p, xo) has the time dependence 

A(p, x ) = A(p)e lp « xa + A(-p)e- ipaX0 (29) 
where po = \p\. Then, the expansion takes the form 

where p ■ x = p^x^ . The transversality condition is satisfied by introducing two 
polarization unit vectors e(p) and 62(39) such that ei • e*2 = ei • P = e*2 • P = 
and ef = ef = 1. Hence if ^4 has to be orthogonal to p, it must be a linear 
combination of ei and e*2, *-e., 

A(P) = J2 ^(p)a a (p) (31) 

a=l,2 

where the factors a Q (p) are complex amplitudes. In terms of a a (p) and a* (p) 
the Hamiltonian looks like a sum of oscillators. 

The passage to the quantum theory is achieved by assigning to each ampli- 
tude a a (p) a Heisenberg operator a(p). Similarly a* (p) maps onto the adjoint 
operator aj^p). The expansion of the vector potential now is 

Ms) = / (2^^ £ e a ( P )[a a (p) e -^- + at(p)e^-] (32) 

with p 2 = and Po = |p|- The operators a a (p) and o) a {p) satisfy commutation 
relations 

[a a (p),al(p')] = 2p (27r) 3 <5(p-p') 
[a a (p),a a ,(p')} = [al(p),al(p')}=0 

(33) 

It is straightforward to check that the vector potential A(x) and the electric 
field E{x) obey the (unconventional) equal-time commutation relation 

[A (x),E f (x')] = -i [s jr ^J^) 5*{x-x') (34) 

where the symbol 1/y 2 represents the inverse of the Laplacian, i.e., the Lapla- 
cian Green function. In the derivation of this relation, the following identity 
was used _ 

£ 4®e£(fl (35) 

a=l,2 P 

These commutation relations are an extension of the canonical ones and it is 
consistent with the transversality condition y • A = 
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In this gauge the (normal-ordered) Hamiltonian is 



6 = I Wjkfo pa £f a(p)&c,m (36) 

The ground state |0) is annihilated by both polarizations a a (p)\0) = 0. The 
single-particle states are dj,(p)|0) and represent photons with momentum p, 
energy pq — \p\ and with the two possible linear polarizations labelled by a = 
1, 2. Circularly polarized photons can be constructed in the usual manner. 



9.3 The Gauge A = 

In this gauge we will apply directly the canonical formalism. In what follows we 
will fix A n = and associate to the three spatial components Aj of the vector 
potential an operator, Aj which acts on a Hilbcrt space of states. Similarly to 
the canonical momentum IX, = —Ej, we assign an operator Tlj. These operators 
obey the equal-time commutation relations 

[A j {x),%,{x')]=i5{x-x')8 jr (37) 

Hence the vector potential A and the electric field E do not commute since 
they are canonically conjugate operators 

[A j (x),E r (x')} = -iS jr 5(x-x') (38) 

Let us now specify the Hilbert space to be the space of states |\E') with wave 
functions which, in the field representation, have the form \&({Aj(x)}). When 
acting on these states, the electric field is the functional differential operator 

sw^im (39) 

In this Hilbert space, the inner product is 

({Aj(x)}\{Aj(x)}}) = Ik;j6 (Aj(x) Aj(x)) (40) 

This Hilbert space is actually much too large. Indeed states with wave functions 
that differ by time-independent gauge transformations 

**({A,-(f)}) = *({A,-(f) - Vi4>W}) (41) 

are physically equivalent since the matrix elements of the electric field operator 
Ej(x) and magnetic field operator Bj(x) — tjkl Vfc Ae(x) are the same for all 
gauge-equivalent states, i.e., 

(^({Mm\E,(m4{Mm) = (niMmEiWmiM*)})) 

(42) 
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The (local) operator Q(x) 

Q(x) = Vj Ej(x) (43) 
commutes locally with the Hamiltonian 

[Q(x),H}=0 (44) 

and, hence, it can be diagonalized simultaneously with H. Let us show now 
that Q(x) generates local infinitesimal time- independent gauge transformations. 
From the canonical commutation relation 

[A j (x),E j ,(x')] = -i5 jj ,S(x-x') (45) 

we get (by differentiation) 

[A j (x),Q(x')} = [A j (x), Vj E j ,(x')]=i V x j S(x-x') (46) 

Hence, we also find 

[i j dz<j>(z)Q(z), Ajix)] = - J dz<t>(z) Vj S(z- x) = Vj<t>($) (47) 



and 



i / dz(f>(z)Q(z) „ —if dz4>(z)Q(z) 
e J ' Aj(x)e J ' = 

-i / dz\7k &(z)E k (z) , i j dz)\j k (j)(z)E k {z) 
= e J Aj(x) e J 

= A 3 (x) + \j 3 (j)(x) 



(48) 



The physical requirement that states that differ by time-independent gauge 
transformations be equivalent to each other leads to the demand that we should 
restrict the Hilbert space to the space of gauge-invariant states. These states, 
which we will denote by |Phys), satisfy 

Q{x) |Phys) = V ' l(f)|Phys) = (49) 

Thus, the constraint means that only the states which obey Gauss' law are in 
the physical Hilbert space. Unlike the quantization in the Coulomb gauge, in 
the A = gauge the commutators are canonical and the states are constrained 
to obey Gauss' law. 

In the Schrodinger picture, the eigenstates of the system obey the Schrodingcr 
equation 



/ 



ds \ [ 'i£w + B > {mnA] = mA] (5o) 



8 



where ^[A] is a shorthand for the wave functional $>({Aj(x)}). In this notation, 
the constraint of Gauss' law is 

V *E 3 (xMA] ee i V J *[A] = (51) 

This constraint can be satisfied by separating the real field Aj (x) into longitu- 
dinal Aj(x) and transverse A J (x) parts 

Aj(x) = Af(x) + Aj(x) = / ^3 ( A f(P) + A J(P)) ( 52 ) 
where Aj(x) and A J (x) satisfy 

VjAj(x) = Af(x) = Vj0(2) (53) 

and cj>(x) is, for the moment, arbitrary. In terms of Af and Aj the constraint 
of Gauss' law simply becomes 

v'jifWf^' (54) 



and the Hamiltonian now is 

6Aj(p)5Aj(-p) SAf(p)SAf(-p) 



*-J*&- 6AZWAZft) ~ SaJsA^p) + ^^)> (55) 



We satisfy the constraint by looking only at gauge-invariant states. Their wave 
functions do not depend on the longitudinal components of A{x). Hence, *£[A] = 
ty[A T ]. When acting on those states, the Hamiltonian is 

Let e*i (p) and e*2 (p) be two vectors which together with the unit vector n p — pj \p\ 
form an orthonormal basis. Let us define the operators (a = 1, 2; j = 1, 2, 3) 



(57) 

These operators satisfy the commutation relations 

[a(p, a),a)(p\ a')} = 6 aa ,5 3 (p - p 1 ) (58) 
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In terms of these operators, the Hamiltonian H and the expansion of the trans- 
verse part of the vector potential are 

/\p\ 
a=l,2 

Aj(x) = / - ^ y €^{p)[a(p,a)e l P- 3 + a\p,a)e-^ 3 ] 
J V(2t) 3 2|p| Q tt 2 



(59) 



We recognize these expressions to be the same ones that we obtained before in 
the Coulomb gauge (except for the normalization factors). 

It is instructive to derive the wave functional for the ground state. The 
ground state |0) is the state annihilated by all the oscillators a(p, a). Hence its 
wave function ^o[A] satisfies 

({A j (x)}\a(p,a)\0)=0 (60) 

This equation is the functional differential equation 

E(P)[^4-=T + \p\Aj(p))M{Aj(p)}) = (61) 
It is easy to check that the unique solution of this equation is 

*o[A] = iVexp[-i j d 3 p\p\Aj(p)Aj(-p)} (62) 
Since the transverse components of Aj(p) satisfy 

A T (rt , t PkA e (p f pxA(p) \ 

A > {f) - eM ^r = \1pY~ ) ( } 

we can write ^o[^4] in the form 

* [A] =AfcM-\J^ (p x MP)) ■ (p x A(-p))} (64) 

It is instructive to write this wave function in position space, i.e., as a functional 
of the configuration of magnetic fields {B(x)}. Clearly, we have 



V x A(p) = -ij -^-^ (vx x A{x)) e 
PXA(-P) = if-^(v x xA(x))e^ 

(65) 
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By substitution of these identities back into the exponent of the wave function, 
we get 

*o[^] = A/"cxp{-iy d 3 x J d 3 x'B(x) ■ B(x')G(x,x')} (66) 
where G(x, x') is given by 

/J3 -ip-(x-x') 
(67 » 

This function has a singular behavior at large values of \p\. We will define a 
smoothed version G\(x, x') to be 

which cuts off the contributions with \p\ ^> A. Also, Ga(x,x') formally goes 
back to G(x,x') as A — > oo . G\(x,x') can be evaluated explicitly to give 

1 f 00 

Gdx,x') = ——— dt sinte-*/^ 5 - 5 



2tt 2 



_1 c^r L 

2 _ ztt\2 1 



(69) 



Thus, 

lim G A (x,f) = - 1 — (70) 

A^oo 2-K Z \X — X'Y 

Hence, the ground state wave functional is 

MA] =AAex P {- i i, / d*x J *x B ^fp ) (71) 
which is only a functional of the configuration of magnetic fields. 



9.4 Path Integral Quantization of Gauge Theories 

We have discuss at length the quantization of the abelian gauge theory (i.e., 
Maxwell's electromagnetism) within canonical quantization in the A = and 
in the Coulomb gauges. Conceptually what we have done is perfectly correct 
although it poses a number of problems. 

1. The canonical formalism is natural in the gauge A — and it can be 
generalized to other gauge theories. However, this gauge is highly non- 
covariant and it is necessary to prove covariance of physical observablcs 
at the end. In addition the gauge field propagator in this gauge is very 
complicated. 
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2. The particle spectrum is most transparent in the transverse (or Coulomb) 
gauge. However, in addition of being non-covariant, it is not possible 
to generalize this gauge to non-Abelian theories due to subtle topologi- 
cal problems (they will be discussed in Physics 583). The propagator is 
equally awful in this gauge. The commutation relations in real space look 
quite different from those in scalar field theory. 

3. In non-Abelian theories, even in the absence of matter fields, the theory 
is already non-linear and needs to be regularized in a manner that gauge 
invariance is preserved. 

4. Although it is possible to use covariant gauges, such as the Lorentz gauge 
d^A^ = 0, the quantization of the theory is these gauges requires a labo- 
rious approach (known as Gupta-Bleuer) of difficult generalization. 

At the root of this problems is the issue of quantizing a theory which has a local 
(or gauge) symmetry in a manner that both Lorentz and gauge invariance are 
kept explicitly. It turns out that path-integral quantization is the most direct 
approach to deal with these problems. 

Let us construct the path integral for the free electromagnetic field. However, 
formally the procedure that we will use can be applied to any gauge theory. We 
will begin with the theory quantized canonically in the gauge A = 0. 

We saw above that, in the gauge Ao = 0, the electric field E is (minus) the 
momentum canonically conjugate to A, the spatial components of the gauge 
field, and obey the equal-time canonical commutation relations 

[Ej{x),A k {x ')} = iS 3 (x-x r ) (72) 

In addition, in this gauge Gauss' Law becomes a constraint on the space of 
states, i.e., 

V • £(aD|Phys) = J (£)|Phys) (73) 

which defines the physical Hilbert space. Here Jo (a;) is a charge density distri- 
bution. In the presence of a set of conserved sources J^(x) (i.e., d^J^ — 0) the 
Hamiltonian of the free field theory is 

H = J d 3 x X - (E 2 + B 2 ) + J d 3 x J ■ A (74) 

We will construct the path-integral in this space. 
Let su denote by Z[J^\ the quantity 

— i I dxoH I —i dxoH \ 
Z[J] = tr'Te J = tr \Te J P (75) 

where tr' means a trace (or sum) over the space of states that satisfy the con- 
straint of Gauss' Law. We implement this constraint by means of the operator 
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P which projects onto these states 



P = Y[s(y-E(x)-J (x)) (76) 

X 

We will now follow the standard construction of the path integral but making 
sure that we only sum over histories that are consistent with the constraint. 
In principle all we need to do is to insert complete sets of states which are 
eigenstates of the field operator A(x) at all intermediate times. These states, 
denoted by \{A(x,x n )}), are not gauge invariant (i.e., they do not satisfy the 
constraint). However, the projection operator P weeds out the unphysical com- 
ponents of these states. Hence, if the projection operator is included in the 
evolution operator, the inserted states actually are gauge- invariant. Thus, to 
insert at every intermediate time Xq (k = 1, . . . , N with N — > oo and Axo — ► 0) 
a complete set of gauge- invariant eigenstates amounts to writing Z[J] as 

N 

Z[J] = H / VAj(x,x k Q ) 
fe=i J 

{{A {x,xl)}\ (l - tAxoH) ]J5 (y • E{x,x k ) - J (x,x k Q )) \{Aj(x , x k +1 )}) 

X 

(77) 

As an operator, the projection operator P is naturally spanned by the eigen- 
states of the electric field operator \{E(x,x )}) 7 i.e., 

Y[ 5 (v • E(x, x ) - Jo(x, so)) = 

X 

j VE(x,x ) \{E(x,x a )})({E(x,x )}\l[S ■ E(x,x ) - J (x,x j) 

X 

(78) 

The delta function has the integral representation 
Y[S (v • E(x,x Q ) - J a (x,x Q )j = 

X 

r iAx / d 3 x A a (x,x a ) (y • E(x,x Q ) - J Q (x, x ) ) 

= Af VA a (x 1 x a )e J v ' 

(79) 
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x e 



Hence, the matrix elements of interest become 

fvAJI ({A(x,x )}\ (l - iAx Q H^j J[S (yjEj - J ) \{A(x,x + Ax )}) 

xo X 

= jvA Q VAVE ]J({A{x,x )}\{E{x,xo)})({E{x,xo)}\{A(x,x + Ax )}) 

f,s , ^ sf-i x T ,~ A ({A(x,xo)}\H\{E(x,x )}) 

/ (Fx Aoix.xQ) (^v • E(x,x ) - Jo(x,x )\ 



iAx 



({A(x,xo)}\{E(x,x )}) 



(80) 



The overlaps are equal to 



i / d 3 x A(x,x a ) ■ E(x,x ) 
: J 



({A(x,x )}\{E(x,x )}) = e 
Hence, we find that the product of the overlaps is given by 

l[ ({A(x,x )}\{E(x,x )}){{E(x,x )}\{A(x,x + Ax )}) = 



(81) 



= e 



dx / d 3 x E(x, x ) ■ doA(x, x ) 



The matrix elements of the Hamiltonian are 
({A(x,x )}\H\{E(x,x )}) 



({A(x,x )}\{E(x,x )}) 



= / d?x 



\(e 2 + Li') + J-A 



1 



(82) 



(83) 



Putting everything together we find that the path integral expression for Z[J] 
has the form 

Z[J] = j VA f tVEe iS [ A K>E] (84) 

where 

VAf, = VAVAq (85) 
and the action S[A^,E] is given by 

r 1 



S[A^,E] 



J d 4 x -E ■ d A - l - (e 1 + B 2 ) - J ■ A + A (y-E- J ) 



(86) 

Notice that the Lagrange multiplier field A , which appeared when we intro- 
duced the integral representation of the delta function, has become the time 
component of the vector potential (that is the reason why I called it ^4 )- 
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Since the action is quadratic in the electric fields, we can integrate them out 
explicitly to find 



/ 



i 

VEe 



J d A x ^E 2 +E- (d A-yA ) 



2 

const, e J 2 K ' ' (87) 

We now collect everything and find that the path integral is 



i J d A x l - (do A- vA, 



/i I d A xC 
VA„ e J (88) 

where the Lagrangian is 

£ = -|f f F""+J/ (89) 

which is what wc should have expected. We should note here that this results is 
valid for all Gauge Theories, Abelian or non-Abelian. In other words, the path 
integral is always the sum over the histories of the field A^ with a weight factor 
which is the exponential of i/h times the action S of the Gauge Theory. 

Therefore we found that, at least formally we can write a functional integral 
which will play the role of the generating functional of the A-point functions of 
this theories, 

-i(0\TA^ 1 (x 1 )...A liN (x N )\0) (90) 



9.5 Path Integrals and Gauge Fixing 

We must emphasize that the expression for the path integral in Eq. (88) is formal 
because we are summing over all histories of the field without restriction. In fact, 
since the action S and the integration measure VA^ are both gauge invariant, 
histories that differ by gauge transformations have the same weight and the 
partition function has an apparent divergence of the form v(G) v , where v(G) 
is the volume of the gauge group G and V is the (infinite) volume of space- 
time. In order to avoid this problem we must implement some sort of gauge 
fixing condition on the sum over histories. We will do so by means of a method 
introduced by L. Faddeev and V. Popov. Although the method works for all 
Gauge Theories, the non-Abelian theories have subtleties and technical issues 
that wc will discus below. We will begin with a general discussion of the method 
and then we will specialize it for the case of Maxwell's theory, the U(l) gauge 
theory without matter fields. 

Let the vector potential A^ be a field which takes values in the algebra of 
a gauge group G, i.e., A^ is a linear combination of the group generators, and 
let U(x) be an unitary-matrix field that takes values on a representation of the 
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group G (please recall our earlier discussion on this subject). For the Abelian 
group U(l), we have 

U(x) = e^ x ) (91) 
where <fi(x) is a real (scalar) field. A gauge transformation is, for a group G 



For the Abelian group U(l) we have 



A^ — Aft + d^c, 



(92) 



(93) 



In order to avoid infinities in Z[J] we must impose restrictions on the sum over 
histories such that histories that are related via a gauge transformation are 
counted exactly once. In order to do that we must find a way to classify the 
vector potentials in classes. We will do this by defining gauge fixing conditions. 
Each class is labelled by a representative configuration and other elements in 
the class are related to it by smooth gauge transformations. Hence, all configu- 
rations in a given class are characterized by a set of gauge invariant data (such 
as field strengths in the Abelian theory) . We must choose gauge conditions such 
that the theory remains local and, if possible, Lorentz covariant. It is essential 
that, whatever gauge condition we use that each class is counted exactly once by 
the gauge condition. It turns out that for the Abelian theory this is trivially the 
case but in non- Abelian theories there are many gauges (such as the Coulomb 
gauge) in which, for topological reasons, a class may be counted more than once. 
(This question is known as the Gribov problem.) Finally we must also keep in 
mind that we are only fixing the local gauge invariance but we should not alter 
the boundary conditions since they represent physical degrees of freedom. 




ations 



Figure 1: The gauge fixing condition selects a manifold of configurations. 
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How do we impose a gauge condition consistently? We will do it in the 
following way. Let us denote by 

g(A^) = (94) 

the gauge condition we wish to impose, where g{A^) is a local diffcrcntiablc 
function of the gauge fields and/or their derivatives, e. g. giA^) = d^A^ 
for the Lorentz gauge or g(A^) = n^A^ for an axial gauge. Note that the 
discussion that follows is valid for all compact Lie groups G of volume v(G). 
For the special case of Maxwell's gauge theory, the gauge group is U(l). Up to 
topological considerations, the group U(l) is isomorphic to the real numbers K, 
even though v(U(l)) = 2ir and v(R) = oo. 

Naively, to impose a gauge condition would mean to restrict the path integral 
by inserting Eq. (94) as a delta function in the integrand, 

Z[J\~ J VA^5{g{A,)) e ^[AJ] (95) 

We will see below that in general this is an inconsistent (and wrong) prescription. 
Following Faddeev and Popov we begin by considering the following integral 

A- 1 ^}^ JvU6(g(AV)) (96) 

where A^ (x) are the configurations of gauge fields related by the gauge trans- 
formation U(x) to the configuration A ll (x), i.e., we move inside one class. 

Let us show that A" 1 ^^] is gauge invariant. We now observe that the 
integration measure VU, usually called the Haar measure, is invariant under 
the composition rule U — ► UU', 

VU = V(UU') (97) 

where U' is and arbitrary but fixed element of G. For the case of G = U(l), 
U = exp(i^) and VU = V<f>. 

Using the invariance of the measure, Eq. (97) we can write 

A" 1 ^'] = JvU5 (g(A u ; u )) = J VU" 6 (<?(<')) = A; 1 ^] (98) 

where we have set U'U = U". Therefore A" 1 [A^] is gauge invariant, i.e., it is 
a function of the dass and not of the configuration A^ itself. Obviously we can 
also write Eq. (96) in the form 

l = A g [A»] JvU6(g(A%)) (99) 

We will now insert the number 1, as given by Eq. (99), in the path integral for 
a general Gauge Theory and find 

Z[J] = JvA^ x 1 x e iS \- A > J } 

= JvA,A g [A,} JvU5(g{A%)) e iS ^ J ^ (100) 
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We now make the change of variables 



A — > A v 



(101) 



where V = V(x) is an arbitrary gauge transformation, and find 

Z[J] =JJvU VAl e^ AV > J ] A g [Al] S (g(A™)) (102) 

(Notice that we have changed the order of integration.) We now choose V = 
?7 _1 , and use the gauge invariance of the action S[A, J], of the measure VA^ 
and of A g [A] to write the partition as 

Z[J] = Jw J VA^ A g [A^} 6 (g{A,)) e iS ^ J ] (103) 

The factor in brackets in Eq. (103) is the infinite constant 

Jt>U = v(G) v (104) 

where v(G) is the volume of the gauge group and V is the (infinite) volume of 
space-time. This infinite constant is nothing but the result of summing over 
gauge-equivalent states. 

Thus, provided the quantity A 9 [^4 M ] is finite and it does not vanish identi- 
cally, we find that the consistent rule for fixing the gauge consists in dividing 
out the (infinite) factor of the volume of the group element but, more impor- 
tantly, to insert together with the constraint S (g(A ll )) the factor A S [A^] in the 
integrand of Z[J], 

Z[J] ~ J VA^ A g [A,] 5 {g(A,)) e iS ^ A (i 05 ) 

We are only left to compute A g [A^] . We will show now that A g [A^\ is a determi- 
nant of a certain operator. The quantity A^L^] is known as the Faddeev-Popov 
determinant. We will only compute first this determinant for the case of the 
Abelian theory U(l). Below we will also discuss the non-Abelian case, relevant 
for Yang-Mills gauge theories. 

We will compute A g [A M ] by using the fact that g[A^] can be regarded as a 
function of U(x) (for A fl (x) fixed). We will now change variables from U to g. 
The price we pay is a Jacobian factor since 



VU = Vg Dct 



6U 



Sg 



(106) 



where the determinant is the Jacobian of the change of variables. Since this is a 
non-linear change of variables, we expect a non-trivial Jacobian. Therefore we 
can write 



A; 1 ^] = J-DU8 (<?«)) = Jvg Det 



SU 



(107) 
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and wc find 



or, conversely 



A^LA^Det 



A g L4„] =Det 



SU 



Sg 
6g 



5U 



9=0 



9=0 



(108) 



(109) 



Thus far all we have done holds for all gauge theories (with a compact gauge 
group). We will specialize our discussion first for the case of the U(l) gauge 
theory, Maxwell's electromagnetism. We will discuss how this applies to non- 
Abelian Yang-Mill gauge theories below. For example, for the particular case of 
the Abelian U(l) gauge theory, the Lorentz gauge condition is obtained by the 
choice giAft) = d^A^. Then, for U(x) = cxp(i(j)(x)) , we get 



giA 1 ^) — (A^ + d^(f>) — d^A^ + d 2 



Hence, 



$g(x) 

5<j>{y) 



d 2 S(x - y) 



(110) 



(111) 



Thus, for the Lorentz gauge of the Abelian theory, the Faddeev-Popov determi- 
nant is given by 

A S L4 M ] =Det9 2 (112) 

which is a constant independent of A^. This is a peculiarity of the Abelian 
theory and, as we will see below, it is not true in the non- Abelian case. 

Let us return momentarily to the general case of Eq. (105), and modify the 
gauge condition from g(A^) = to g{A^) = c(x) 1 where c(x) is some arbitrary 
function of x. The partition function now reads 

Z[J] ~ J VA„ A g [A^} 6 {g(A,) - c{x)) e lS ^ J ] (113) 

We will now average over the arbitrary functions with a Gaussian weight (prop- 
erly normalized to unity) 



Z a [J] = Af J VAf, Vce 1 1 



d±x C(X)2 



2" A g [A,}S(g(A,)-c(x)) e iS \ A ^ 



Af / VA„ A g (A fl ) e 



+i J d A x 



C[A,J]-^-(g(A,)) 2 



(114) 



From now on we will restrict our discussion to the U(l) Abelian gauge theory 
(the electromagnetic field) and giA^) = d^A^. From Eq. (114) we find that in 
this gauge the Lagrangian is 



C a = -\Fl v -J^-±-(d,A») 2 



(115) 
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The parameter a labels a family of gauge fixing conditions known as the Feynman- 
't Hooft gauges. For a — > we recover the strong constraint d^A^ = 0, the 
Lorentz gauge. From the point of view of doing calculations the simplest is the 
gauge a — 1 (the Feynman gauge) as we will see now. After some algebra is 
straightforward to see that, up to surface terms, the Lagrangian is equal to 



r — - a 

*~a — 2 v 



a 



A-v J ^A^ 



(116) 



and the partition function reduces to 



Z[J] = M Dot [d 2 ] J VA 
Hence, in a general gauge labelled by a, we get 



J d 4 x£ a [A,J] 



(117) 



Z[J] = N Dct [d 2 ] Det 



t v d 2 



a-l 



a 



-1/2 



x e 



l -j 'd 4 xj d 4 y J„{x) G^(x - y) J v (y) 



(118) 



where G^(x — y) is the propagator in that gauge. 

The form of Eq. (118) may seem to imply that Z[J] is gauge dependent. This 
cannot be correct since the path integral is by construction gauge-invariant. We 
will show in the next subsection that gauge invariance is indeed protected. This 
result comes about because J M is a conserved current, and as such it satisfies 
the continuity equation <9 M J M = 0. For this family of gauges, the propagator 
takes the form 



G^(x - y) 



a - 1 d»d v 



a 



d 2 



G(x - y) 



(119) 



where G(x — y) is the propagator of the scalar field. 

Thus, as expected for a free field theory, Z[J] is a product of two factors: 
a functional (or fluctuation) determinant, and a factor that depends solely on 
the sources J M which contains all the information on the correlation functions. 
For the case of a single scalar field we also found a contribution in the form of a 
determinant factor but its power was —1/2. Here there are two such factors. The 
first one is the Faddeev-Popov determinant. The second one is the determinant 
of the fluctuation operator for the gauge field. However, in the Feynman gauge, 
a = 1, this operator is just g^d 2 , and its determinant has the same form as 
the Faddeev-Popov determinant except that it has a power —4/2. This is what 
one would have expected for a theory with four independent fields (one for each 
component of A^). The Faddeev-Popov determinant has power +1. Thus the 
total power is just 1 — 4/2 = —1, which is the correct answer for a theory with 
only two independent fields. 
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9.6 The Propagator 

For general a, G fJ/l/ (x — y) is the solution of the Green's function equation 



a 



G vX {x -y) = g%5 {x - y) 



(120) 



Notice that in the special case of the Feynman gauge, a = 1, this equation 
becomes 

d 2 G^(x-y) = g^5 4 (x-y) 
Hence, in the Feynman gauge, G^ v (x — y) takes the form 



C? v (x-v)=g» v G(x-y) 
where G(x — y) is just the propagator of a free scalar field, i.e., 

8 2 G(x -y)= 5\x - y) 
However, in a general gauge the propagator 

G^(x -y) = -i{Q\TA^x)A u {y)\0) 



(121) 
(122) 

(123) 

(124) 



does not coincide with the propagator of a scalar field. Therefore, G^ v {x — y), 
as expected, is a gauge dependent quantity. In spite of that it does contain 
physical information. Let us examine this issue by calculating the propagator 
in a general gauge a. 

The Fourier transform of G l _ lv (x — y) in D space-time dimensions is 



This a solution of Eq. (120) provided G^p) satisfies 

G»\(p) = 3a 



-g^p 2 + - — -p"p" 



a 



The formal solution is 



G>„(p) 



P 



a — 1 pVp v 
a p 2 



(125) 



(126) 



(127) 



In space-time the form of this (still formal ) solution is given by Eq. (119). 
In particular, in the Feynman gauge a = 1, we get 



whereas in the Lorentz gauge we find instead 



y ^2 



(128) 



(129) 
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Hence, in all cases there is a pole in p 2 in front of the propagator and a matrix 
structure that depends on the gauge choice. Notice that the matrix in brackets 
in the Lorentz gauge, known as the transverse projection operator, satisfies 



p2 



(130) 



which follows from the gauge condition d^A^ = 0. 

The physical information of this propagator is condensed in its analytic 
structure. It has a pole at p 2 = which implies that p a = = \p\ is the 

singularity of G ' M „ (p) . Hence the pole in the propagator tells us that this theory 
has a massless particle, the photon. 

To actually compute the propagator in space-time from G flu (p) requires that 
we define the integrals in momentum space carefully. As it stands, the Fourier 
integral Eq. (125) is ill defined due to the pole in G^ v (p) at p 2 = 0. A proper 
definition requires that we move the pole into the complex plane by shifting 
p 2 — > p 2 + ie, where e is real and e — > + . This prescription yields the Feynman 
propagator. We will see in the next section that this rule applies to any theory 
and that it always yields the vacuum expectation value of the time ordered 
product of fields. For the rest of this section we will use the propagator in the 
Feynman gauge which reduces to the propagator of a scalar field. This is a 
quantity we know quite well, both in Euclidean and Minkowski space-times. 

9.7 Physical meaning of Z[J\ and the Wilson loop operator 

We discussed before that a general property of the path integral of any theory 
is that, in imaginary time, Z[0] is just 

Z[0] = (0|0) ~ e~ TE ° (131) 

where T is the time span (i.e., T — > oo) (watch out, here T is not the tem- 
perature!), and Eq is the vacuum energy Thus, if the sources J M are static (or 
quasi-static) we get instead 

m^ e -T[E (J)-E ] (132) 

Thus, the change in the vacuum energy due to the presence of the sources is 

U(J) = E (J) -E = T lim In ||^ (133) 

As we will see, the behavior of this quantity has a lot of information about the 
physical properties of the vacuum (i.e., the ground state) of a theory. Quite 
generally, if the quasi-static sources J M are well separated from each other, U (J) 
can be split into two terms: a self-energy of the sources, and an interaction 
energy, i.e., 

U(J) = E seli -energy [J]+V int [J] (134) 
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As an example, we will now compute the expectation value of the Wilson 
loop operator, 



W r = (0\Te Jr 



dx^ 



|0> 



(135) 



where T is the closed path in space-time shown in the figure. Physically, what 
we are doing is looking at the electromagnetic field created by the current 



(136) 



where is the set of points of space-time on the loop T , and is a unit 
vector field tangent to T. The loop T has time span T and spatial width R. 
We will be interested in loops such that T >• R so that the sources are turned 
on adiabatically in the remote past and switched off also adiabatically in the 
remote future. By current conservation the loop must be oriented. Thus, at 
a fixed time xq the loop looks like a pair of static sources with charges ±e at 
±R/2. In other words we are looking of the affects of a particle-antiparticle 
pair which is created at rest in the remote past, the members of the pair are 
slowly separated (to avoid bremsstrahlung radiation) and live happily apart 
from each other, at a prudent distance R, for a long time T, and are finally 
(and adiabatically) annihilated in the remote future. Thus, we are in the quasi- 
static regime described above and Z[J]/Z[0] should tell us what is the effective 
interaction between this pair of sources ( "electrodes" ) . 



adiabatic 
switching off 




Figure 2: The Wilson loop operator can be viewed as representing a pair of 
quasi-static sources of charge ±e separated a distance R from each other. 
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What arc the possible behaviors of the Wilson loop operator in general 
(that is for any gauge theory)? The answer to this depends on the nature of 
the vacuum state. We will se that a given theory may have different vacua 
or phases (as in thermodynamic phases) and that the behavior of the physical 
observables is different in different vacua (or phases) . Here we will do an explicit 
computation for the case of Maxwell's U(l) gauge theory. However the behavior 
that will find only holds for a free field and it is not generic. What are the 
possible behaviors, then? A loop is an extended object (as opposed to a local 
field operator) characterized by its geometric properties: its area, perimeter, 
aspect ratio, and so on. We will show later on that these geometric properties 
of the loop to characterize the behavior of the Wilson loop operator. Here are 
the generic cases: 

1. Area Law: Let A = RT be the area of the loop. One possible behavior of 
the Wilson loop operator is the area law 

W r ~ e~ aRT (137) 

w We will show later on that this is the fastest possible decay of the Wilson 
loop operator as a function of size. The quantity a is known as the string 
tension. If the area law is obeyed the effective potential for R large (but 
still small compared to T) behaves as 

V int {R) = lim -llnVKr = <tR (138) 

T — >oo I 

Hence, in this case the energy to separate a pair of sources grows linearly 
with distance and the sources are confined. We will say that in this case 
the the theory is in confined phase. 

2. Perimeter Law: Another possible decay behavior (weaker than the area 
law) is a perimeter law 

W r ~e-P( R + T )+0(e- R /t) (139) 

where p is a constant with units of energy, and £ is a length scale. This 
decay law implies that in this case 

Mnt ~ P + const. e- R /Z (140) 

Thus, in this case the energy to separate two sources is finite. This is a 
deconfined phase. However since it is massive (with a mass scale m ~ £ _1 ) 
there are no long range gauge bosons. This is the Higgs phase. It bears a 
close analogy with a superconductor. 

3. Scale Invariant: Yet another possibility is that the Wilson loop behavior 
is determined by the aspect ratio R/T or T/R, 

-"(- + -) 
W T ~e \ T R ) (141) 
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where a is a dimensionless constant. This behavior leads to an interaction 

Vint ~ ~ (142) 

which coincides with the Coulomb law in 4 dimensions. We will see that 
this is a deconfined phase with massless gauge bosons (photons). 

We will now compute the expectation value of the Wilson loop operator in 
Maxwell's U(l) gauge theory. We will return to the general problem when we 
discuss the strong coupling behavior of gauge theories. We begin by using the 
analytic continuation of Eq. (118) to imaginary time, i.e., 

m = J ^ OT -.^sW* , ' w ^ (I) * M> * ( ' , 

= H^[ 3 r i e-^l' iX 'i dyAA ' {X)AM) (143) 

where (A IM (x)A v (y)) is the Euclidean propagator of the gauge fields in the family 
of gauges labelled by a. In the Feynman gauge a = 1 the propagator is given 
by the expression 

(M*)Mv)) = <W / (0j £ ' {x » ^ (144) 

where \i = 1, . . . , D. After doing the integral we find that the Euclidean propa- 
gator (the correlation function) in the Feynman gauge is 



F 



D 



(A^(x)A v (y)) = V ^o/a \ x _ Jd-2 (14:,) 
Therefore E[J] — Eo is equal to 

' D -l 



e 2 f f ,„ V 2 



E[J] — Eq = lim — f f dx ■ dx 



r^ooZr/r/r A-k d ' 2 \x - y\ D ~ 2 

r( D 1 

g 2 r+T/2 r+T/2 \~2 



= 2 x self - energy - — 2 / dx D I dy D 

Z1 J-T/2 J-T/2 



,D-2 



T/2 J-T/2 4*7T ' \x — y\ 

(146) 
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where \x — y\ 2 = (xd — Vd) 2 + R 2 - The integral in Eq. (146) is equal to 



+T/2 r+T/2 

dx D / dy D 



D-2 



T/2 J-T/2 47T ' |x — y 

+T/2 , + (T/2-,)/fl ^ 1 T ,D-2S 

ds / — 77- oW „ „ n 5 X 



T/2 J-(T/2+,)/R (^ + 1)^-2)72^-3 4^/2 

+T/2 _ ^+00 ^ p (Z2_2) 



1 /-m/^ r- 

— It/2 ds 



^ D ~ 3 J-t/2 J-co (t 2 + iy D -v/ 2 4tt B / 2 

r(^) 47r C / 2 



(147) 



where 



I» = / dtf^e-t (148) 
Jo 

is the Euler Gamma function. In Eq. (147) we already took the limit T/R — > oo. 
Putting it all together we find that the interaction energy of a pair of static 
sources of charges ±e separated a distance R in D dimensional space-time is 
given by 

r(^=i) e 2 

Mnt(i?) = " 2^ d -^)/\d - 3) (U9) 

This is the Coulomb potential in D space-time dimensions. In particular, in 
D = 4 dimensions we find 

v^R) = - (£) ± d50) 

where the quantity e 2 /47r is the fine structure constant. 

Therefore we find that, even at the quantum level, the effective interaction 
between a pair of static sources is the Coulomb interaction. This is true because 
Maxwell's theory is a free field theory. It is also true in Quantum Electrodynam- 
ics (QED), the Quantum Field Theory of electrons and photons, at distances R 
much greater than the Compton wavelength of the electron. However it is not 
true at short distances where the effective charge is screened by fluctuations of 
the Dirac field and the potential becomes exponentially suppressed. In contrast, 
in Quantum Chromodynamics (QCD) the situation is quite different: even in 
the absence of matter, for R large compared with a scale £ determined by the 
dynamics, the effective potential grows linearly with R. This effect in known as 
confinement, and the scale £ is the confinement scale. Conversely, the potential 
is Coulomb-like at short distances, a behavior known as asymptotic freedom. 



26 



9.8 Path-Integral Quantization of Non- Abelian Gauge The- 
ories 

Most of what we did for the Abelian case, carries over to non-Abelian gauge 
theories where, as we will see, it plays a much more central role. In this section 
we will discuss the general properties of the path-integral quantization of non- 
Abelian gauge theories, but we will not deal with the non-linearities here. 

The path integral Z[J] for a non-Abelian gauge field = A^\ a in the 
algebra of a simply connected compact Lie group G, whose generators are the 
Hermitian matrices A a , with gauge condition(s) g a [A] is 

Z[J] = J VAl e iS ^ J] S (g[A}) App [A] (151) 

where we will use the family of covariant gauge conditions 

g a [A] = d^A^x) + c a (x) = (152) 

and where App[A] is the Faddeev-Popov determinant. Notice that we impose 
one gauge condition for each direction in the algebra of the gauge group G. 
We will proceed as we did in the Abelian case and consider an average over 
gauges. In other words we will work in the manifestly covariant Feynman-'t 
Hooft gauges. 

Let us work out the structure of the Faddeev-Popov determinant for a general 
gauge fixing condition g a L4]. Let U be an infinitesimal gauge transformation, 

U~ l + ie a (x)X a + ... (153) 

Under a gauge transformation the vector field A^ becomes 

A u ^=UA^U- 1 + i{d^U)U- 1 = A^ + 8A^ (154) 

For an infinitesimal transformation the change of A^ is 

SAp = ie a [A a , A„] - d^ a \ a + 0(e 2 ) (155) 

where A a are the generators of the algebra of the gauge group G. 
On the other hand, since 

5g a dg a 5A C 



M 

5e b dA c u Se b 



(156) 



we can also write 



8A<^ = 2ie b tr (A c [A fc , A^)-2 d^ b tr (A c A b ) + 0{e 2 ) 

= ie\i(\ c [A fc ,A d ])^-d M e b <5 6c 

= -2 f bde e"tr (A c A e ) A d - d^e b S bc 

= -f bde e» A* - d^e» S bc (157) 
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where f abc are the structure constants of the Lie group G. 
Hence, we find 

8A C (x) 

-J^y = [d,S bc + f bcd Af] 6(x -y) = -B?[A]6{x y) (158) 

where we have denoted by [A] the covariant derivative in the adjoint repre- 
sentation, which in components is given by 

Df[A] = 5 ab dp - f abc A^ (159) 

Using these results we can put the Faddeev-Popov determinant (or Jacobian) 
in the form 

S A ^Jdg a 5A^ 



A^]=Det^j=Det^_^ (160 ) 

We will now define an operator Mpp whose matrix elements are 

dg SA C 

(x,a\M FP \y,b) = (x,a\-^-^\y,b) 
dg a (x) SA^z) 



dAfa) Se b (y) 



For the case of g a [A] = d^A"(x) — c a (c), appropriate for the Feynman-'t Hooft 
gauges, we have 

^ -«-•"(*-*) < 162 » 

and also 

(x,a\M FP \y,b) = - [ 6 ac d%6(x - z)Df[A]6(z - y) 

J z 

8 ac S(x - z) d?D cb S(x - y) 

= -d^D ab 5{x - y) (163) 

Thus, the Faddeev-Popov determinant now is 

App = Det (d^D^A]) (164) 

Notice that in the non-Abelian case this determinant is an explicit function of 
the gauge field. Since it is a determinant, it can be written as a path integral 
over a set of fcrmionic ghost fields, denoted by r) a (x) and fj a (x), one per gauge 
condition (i.e. one per generator): 



r i f d D xfj a (x)d^Df[A]r lb (x) 

Det = J VvaVfja e J 



(165) 
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Notice that these are fields quantized with the "wrong" statistics. In other 
words, these "particles" do not satisfy the general conditions for causality and 
unitarity to be obeyed. Hence these ghosts cannot create physical states (thereby 
their ghostly character). 

The full form of the path integral of a Yang-Mills gauge theory with coupling 
constant g, in the Feynman-'t Hooft covariant gauges with gauge parameter A, 
is given by 

r i / d D x CymIAViV] 

Z = VAVr]Vfj e J (166) 

where Cym is the effective Lagrangian density 

C YM [A, V ,fj\ = --LtrF^F"" + A {d^f - fj 8^[A] r, (167) 

Thus the pure gauge theory, even in the absence of matter fields, is non-linear. 
We will return to this problem later on when we look at both the perturbative 
and non-perturbative aspects of Yang-Mills gauge theories. 



10 BRST Invariance 

In the previous section we developed in detail the path-integral quantization of 
non-Abelian Yang-Mills gauge theories. We payed close attention to the role 
of gauge invariance and how to consistently fix the gauge in order to define 
the path-integral. Here we will show that the effective Lagrangian of a Yang- 
Mills gauge field, Eq. (167), has an extended symmetry, closely related to 
supersymmetry. This extended symmetry plays a crucial role in proving the 
renormalizability of non-Abelian gauge theories. 

Let us consider the QCD Lagrangian in the Feynman-'t Hooft covariant 
gauges (with gauge parameter A and coupling constant g). The Lagrangian 
density C of this theory is 

C = ^{ip-rn}^- \f; v FZ v - ^B a B a + B a &>Al - f} a d»D a J> v b 

(168) 

Here ip is a Dirac Fermi field, which represents quarks and it transforms under 
the fundamental representation of the gauge group G; the "Hubbard-Stratonovich" 
field B a is an auxiliary field which has no dynamics of its own and it transforms 
as a vector in the adjoint representation of G. 

Becchi, Rouet, Stora and Tyutin realized that this gauge-fixed Lagrangian 
has the following ( "BRST" ) symmetry, where e is an infinitesimal anti-commuting 
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parameter: 





= eD; b r, b 


(169) 




= iger) a t a ip 


(170) 


<5r/ a 


= -\ gef abc VbVc 


(171) 


5fj a 


= eB a 


(172) 


SB a 


= 


(173) 



Here Eq. (169) and Eq. (170) are local gauge transformations and as such leave 
invariant the first two terms of the effective Lagrangian C of Eq. (168). The 
third term of Eq. (168) is trivial. The invariance of the fourth and fifth terms 
holds because the change of 5A in the fourth term cancels against the change 
of f\ in the fifth term. Finally, it remains to see that the changes of the fields 
and r\ in the fifth term of Eq. (168) cancel out. To see that this is the case 
we check that 



S {Dfr?) = Dp^+gf^SA^ = -^ghf^f^ {A^rf + A^rf + A^rf) 

(174) 

which vanishes due to the Jacobi identity for the structure constants 



jade jbed _|_ jbde jcad _j_ jede jabd q (175) 

or, equivalcntly from the nested commutators of the generators t a : 

[t°, [t b ,t c ]] + [t h ,[t c ,t a \] + [t c , [t a ,t b ]] =0 (176) 

Hence, BRST is at least a global symmetry of the gauge-fixed action with gauge 
fixing parameter A. 

This symmetry has a remarkable property which follows from its fermionic 
nature. Let (ft be any of the fields of the Lagrangian and Q<fi be the BRST 
transformation of the field, 

6(f> = eQ<f> (177) 

For instance, 

Q a Al = Dfrf (178) 
and so on. It follows that for any field </> 

Q 2 cf) = (179) 



i.e. the BRST transformation of Q(j> vanishes. This rule works for the field A^ 



due to the transformation property of S(D° l b r] b ). It also holds for the ghosts 



since 

QV = ^g 2 f abc .f bde V c V d V e = (180) 



which holds due to the Jacobi identity. 
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What are the implications of the existence of BRST as a continuous sym- 
metry? To begin with it implies that there is a conserved self-adjoint charge 
Q that must necessarily commute with the Hamiltonian H of the Yang-Mills 
gauge theory. Above we saw how Q acts on the fields, Q 2 (j> = 0, for all the fields 
in the Lagrangian. Hence, as an operator Q 2 = 0, that is, the BRST charge 
Q is nilpotent, and it commutes with H. Let us now show that Q divides the 
Hilbert space of the eigenstates of H is three sectors 

1. Many eigenstates of H must be annihilated by Q for Q 2 = to hold. Let 
Hi be the set of eigenstates of H which are not annihilated by Q. Hence, 
if |V>i) € Hi, then Q\ipi) ^ 0. Thus, the states in Hi are not BRST 
invariant. 

2. Let us consider the subspace of states H2 of the form ^2) = Q\ipi), 

H2 = QHi- Then, for these states Q\ip2) = Q 2 \^i) = 0. Hence, the states 
in H2 arc BRST invariant but are the BRST transform of states in Hi- 

3. Finally, let Ho be the set of eigenstates of H that are annihilated by Q, 
Q\ipo) = 0, but which are not in H2, i-e. \ipo) 7^ Q\ipi}- Hence, the states 
in Ho arc BRST invariant and are not the BRST transform of any other 
state. This is the physical space of states. 

It follows from the above classification that the inner product of any pair of 
states in H2, ^2) and |^> 2 )> have zero inner product: 

W*|V>2> = <^l|W2>=0 (181) 

where we used that ^2) is the BRST transform of a state in Hi, \ipi). Similarly, 
one can show that if \tpo) € Ho, then (ip2\i>o) = 0. 

What is the physical meaning of BRST and of this classification? Peskin 
and Schroeder give a simple argument. Consider the weak coupling limit of 
the theory, g — > 0. In this limit we can find out what BRST does by looking 
at the transformation properties of the fields that appear in the Lagrangian of 
Eq. (168). In particular, Q transforms a forward polarized (i.e. longitudinal) 
component of into a ghost. At g — 0, we see that Qn = and that the anti- 
ghost fj transforms into the auxiliary field B. Also, at the classical level, B = 
Ad^Afi. Hence, the auxiliary fields B are backward (longitudinally) polarized 
quanta of A^. Thus, forward polarized gauge bosons and anti-ghosts arc in 
Hi, since they are not the BRST transform of states created by other fields. 
Ghosts and backward polarized gauge bosons are in H2 since they are the BRST 
transform of the former. Finally, transverse gauge bosons are in Ho- Hence, 
in general, states with ghosts, anti-ghosts, and gauge bosons with unphysical 
polarization belong either to Hi or H2- Only the physical states belong to Ho- 
lt turns out that the S-matrix, when restricted to the physical space Ho, is 
unitary (as it should). 
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